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Abstract 

The "body fixed frame" witli respect to local gauge transformations is intro- 
duced. Rigid gauge "rotations" in QCD and their Schrodinger equation are studied 
for static and dynamic quarks. Possible choices of the rigid gauge field configuration 
corresponding to a nonvanishing static colormagnetic field in the "body fixed" frame 
are discussed. A gauge invariant variational equation is derived in this frame. For 
large number N of colors the rigid gauge field configuration is regarded as random 
with maximally random probability distribution under constraints on macroscopic- 
like quantities. For the uniform magnetic field the joint probability distribution of 
the field components is determined by maximizing the appropriate entropy under the 
area law constraint for the Wilson loop. In the quark sector the gauge invariance 
requires the rigid gauge field configuration to appear not only as a background but 
also as inducing an instantaneous quark-quark interaction. Both are random in the 
large N limit. 



1 Introduction 

Studies of non perturbative aspects of dynamics of non abelian gauge fields will con- 
tinue to remain one of the focuses of theoretical activities. These fields appear at all 
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levels of the "elementary" interactions and even begin to enter at a more plienomeno- 
logical macroscopic level in condensed matter systems. Quantum Cliromodynamics 
represents a prime example of a strongly coupled theory with non abelian gauge fields. 
Despite many efforts, e.g. instantons large N expansion [0, H, ^, ^, ||, lattice 
gauge theory and strong coupling expansion 0, |], topological considerations |10|, p!T|| , 
QCD sum rules [0, "spaghetti" vacuum |]13|, light cone approach ||14|, explicit color 
projection and others [jl6|, 0, |18|, the quantitative understanding of the basic 



QCD features is still far from satisfactory. A sustained effort with different angles of 
attack is clearly in order with the hope that accumulated qualitative experience will 
finally lead to the development of quantitative calculational tools. This paper is a 
contribution to this effort. 

Invariance under local gauge transformations is the most important feature of a 
non abelian gauge theory. In the framework of the Hamiltonian formulation of QCD 
I wish to explore the consequences of this invariance using some general methods 
common to molecular and nuclear physics. I wish to define an appropriate general- 
ization of the body-fixed (intrinsic, rotating) frame formalism in the context of the 
local gauge transformations. After doing so one can attempt to separately investigate 
the dynamics of the gauge "rotations" of the frame and the intrinsic frame dynamics. 
These would be the analogs of rigid rotations and intrinsic vibrations in molecular 
and nuclear physics. Most of my interest in this paper will concentrate on the study 
of the "rigid gauge rotations". The study of the couplings between the "rotations" 
and the "vibrations" of the gauge field is deferred to future work. To avoid misunder- 
standing I wish to stress that by "gauge rotations" or "rotations of the gauge field" 
in this paper I will always mean Eq. below, which includes the proper SU(N) 
rotation as well as the inhomogeneous "shift" term. 

Perhaps the most important conceptual advantage of using the body-fixed frame 
associated with a given symmetry lies in the fact that one can freely approximate 
the dynamics in this frame without fears to violate the symmetry. In particular the 
use of this formalism appears to be fruitful provided there exist such a body-fixed, 
intrinsic frame in which the "rotational - vibrational" coupling can be considered as 
small. This in turn generically happens when the "rotational inertia" is much larger 
than the "inertia" associated with the intrinsic motion so that a variant of the Born- 
Oppenheimer approximation is valid. A typical situation is when the system's ground 
state is strongly "deformed" away from a symmetric state. By deformation I mean 
absence of symmetry with respect to "intrinsic transformations", i.e. transformations 
in the body- fixed frame and not with respect to the transformations in "laboratory". 
Absence of symmetry in "laboratory" would correspond to the symmetry breakdown 
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which can not occur for a local gauge symmetry. Quantum mechanical examples of 
deformed bodies are e.g. non spherical molecules, deformed nuclei, etc. 

I do not have a priori arguments that the QCD vacuum is strongly "deformed" 
in the above sense. Appearance of various QCD condensates, Ref. [0, suggests 



that this may be true. The condensate wavefunction should then play a role of the 
strongly deformed configuration. Another positive indication is the large N master 
field concept, Ref. 0, according to which a special gauge field configuration should 
exist which dominates the vacuum wave function or the corresponding functional 
integral. It is expected, however, that the master field is not simply a fixed classical 
configuration. It should rather be regarded as a statistical distribution allowing to 
calculate quantities which are analogous to macroscopic thermodynamic quantities in 
statistical physics, i.e. such that their fluctuations are suppressed in the large N limit, 
Refs. [p!9| , pO|] . Glimpses of the meaning of these vague notions were found in various 
matrix models, cf. Refs. ^, and, e.g. in 1+1 dimensional QCD, Ref. PT[] . 

If this view point is correct then suitably chosen rigidly rotating "deformed" gauge 
field could play a role of the master field provided one understands in which sense it 
should also be statistical. The following formalism will clarify some of these issues 
and provide a general framework in which they could be further discussed. 

Works in the spirit of our study have already appeared in the past,cf. Refs. 



2^, |2^, ^ and the analogy with various types of rotational motion is frequently 
used in QCD. In this sense the present study is a continuation of these works. 

This paper is organized as following. In Section 2 I introduce the transformation 
to the body-fixed frame in the context of the simplest model of the gauge rotational 
motion - the rigid gauge rotor. Giving a natural definition of the rigid gauge "rota- 
tions" I proceed to determine the appropriate generalization of the standard space 
rigid rotor results - the moment of inertia tensor, the "body-fixed" frame, the genera- 
tors of the "body- fixed" gauge group in terms of which the character of "deformation" 
can be classified, etc. Despite severe limitation on the set of the allowed gauge field 
configurations the model is gauge invariant. I work out in Section 3 the quantum 
mechanics of the model. As with the space rotor the generators of the "laboratory" 
and "body- fixed" gauge transformations provide a complete set of quantum numbers 
for the wave functionals of the model. The vacuum has zero energy and is the most 
disordered state. Higher states correspond to the presence of very heavy, i.e. static 
quarks and antiquarks in the system. As an important example I consider the wave 
function and the corresponding Schrodinger equation for a pair of static quark and 
antiquark. This and other similar equations in the model are simple matrix equation 
with the inverse "moment of inertia" determining the interaction between the color 
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sources and depending on the assumed rigid gauge configuration which plays the role 
of the " free parameter" . 

In Section 4 I discuss the meaning of the results obtained so far and possible 
choices of the rigid gauge field which physically represents a non vanishing colormag- 
netic field in the body- fixed frame. In the ground state this frame does not "rotate" 
but has random orientations in local color spaces at every space point. Introduction 
of static quarks forces the frame to "rotate" quantum mechanically at the points 
where the quarks are situated. The energy eigenvalues of these "rotations" are the 
energies of the quantum states of the colorelectric field generated by the quarks. The 
propagator of this field is the moment of inertia of the model and depends explicitly 
on the assumed configuration of the rigid static colormagnetic field. For zero field the 
propagator is a simple Coulomb while for a uniform field diagonal in color the prop- 
agator behaves asymptotically as a decaying Gaussian. The so called dual Meissner 
effect picture of the confinement, Ref. could be implemented if a configuration 
of the rigid colormagnetic field is found which "channels" the colorelectric field and 
makes its propagator effectively one dimensional. It turns out that the creation of 
such a magnetic "wave-guide" is connected with existance of a zero eigenvalue of a 
certain operator in the model. 

Since the quark color degrees of freedom are treated quantum mechanically the 
model allows for a possibility that confinement of fundamental representations does 
not automatically mean confinement of higher representations. I discuss this possi- 
bility and derive a variational equation for the rigid field. This equation is fully gauge 
invariant. 

In Section 5 expecting that rigid gauge rotations should be relevant for the master 
field concept I study the model in the large N limit. Any candidate for the master 
field must be allowed to undergo free "gauge rotations" which can not be frozen by 
this limit and should induce an interaction between the quarks. Going to the body 
fixed frame of these rotations I regard the rigid gauge field configuration as random 
and introduce a natural requirement that it is least biased under constraints that it 
should reproduce gauge invariant quantities which can be regarded macroscopic-like 
in the large N limit. This means that it should be maximally random under these 
constraints. In order to make these ideas explicit I discuss in some detail the case of 
the uniform colormagnetic field. Such configuration in QCD was already discussed in 



the past, Refs. |13[ but it seems that its appearance in the interaction is a novel 



feature of the model. The detailed form of this interaction depends on the differences 
of the color components of the magnetic field. It is not confining for any finite number 
of colors. For N —>■ ool assume that the form of the density of the color components of 
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the field is known. In 2+1 dimensions I choose it such that it gives area law for space 
oriented Wilson loops. I treat then the entire distribution of these components as a 
joint distribution of their probabilities and regard the adopted "single component" 
density as an analog of a macroscopic quantity that must be reproduced by this joint 
distribution. I postulate that it must otherwise be maximally random, i.e. must have 
the maximum entropy (minimum information content) under suitable constraints. 
In this way I derive the maximally random distribution for this model. I discuss its 
relation to the large N limit of the Schrodinger equation for the static quark-antiquark 
system. I also give possible generalizations of this development to 3+1 dimensions. 

In Section 6 I include dynamical quarks and show that the rigid gauge rotor limit 
corresponds exactly to the limit in QED in which only the instantaneous Coulomb 
interaction between the charges is retained. The major difference in QCD is that 
in this limit the quarks not only interact instantaneously via a more complicated 
interaction controlled by the rigid gauge field configuration, but at the same time are 
also found in a static colormagnetic field induced by this configuration. This dual 
appearance of the rigid field is a consequence of the gauge invariance and in the large 
N limit is apparently the way the master field should enter the quark sector of the 
theory. According to the ideology developed in Section 5 both the field in which the 
quarks move and their interaction should be considered as random in the large N 
limit. The random interaction between the quarks opens interesting possibilities to 
discuss the relationship between confinement and localization. 

The body fixed Hamiltonian with dynamical quarks is gauge invariant. Its invari- 
ance with respect to global symmetries however is not guaranteed for an arbitrary 
choice of the rigid gauge configuration. I discuss possible variational approaches to 
determine this configuration and derive an analogue of the Hartree-Fock equations 
for the model. 

In the rest of the Introduction I will establish my notations, cf., Ref. 
the QCD Hamiltonian in d=3 space dimensions in the = gauge, 

H = y d^x[{El{x)r + (Biix))'] + J d'xq;{x)[a^ (^p^ - gK{x)^^ + /5m]g,(x). 

(1.1) 

with 

Bl{x)=e,,u{d,Al + gU,J^,A',), (1.2) 

i,j,k = 1,...,3; 7,(5=1,...,N; a = 1,...,N'^ - 1 for SU(N) gauge group and fate - 
the structure constants of the SU(N). Dirac and flavor indices are omitted and the 
summation convention for all repeated indices is employed here and in the following. 



2£]. I consider 
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The gluon vector potential Al^{x) and minus the electric field — -E* (a;) are canonically 
conjugate variables, 

[Eiix),Aiiy)]=t5,,5,,5{x-y), (1.3) 

and the quark fields obey the standard anticommutation relations. 

The Hamiltonian ( |1 . 1|) is invariant under the time independent gauge transforma- 
tions. Using the matrix valued hermitian fields 

A'^pi^) = Kix)^, Kpi^) = Kix)^, (1-4) 

where A" are the SU(N) generators with the properties 

[A^A''] = 2zUcX''; {A^A^} = ^(5,6 + 24fecA=; 
Tr A-^A^ = 25,,; X-^X-^ = 2[6^sSp,-j^6^p6,s] (1.5) 
one can write the gauge transformation as 

A' ^ SA'S+ + -Sd'S+; E'^SE'S+; q ^ Sq (1.6) 
9 

where S(x) are time independent but x - dependent unitary N xN matrices, elements 
of the SU(N) group. The generators of this transformation 

G,{y)^Gt{y) + Gl{y), (1.7) 

G^{y) = d,EM + gU,Al{y)Eliy), Gl{y) = -gq^{y)^q{y) (1.8) 
are conserved, 

^^ = ^[H,Gaix)] = (1.9) 
and it is consistent to impose the Gauss law constraints 

G,(x)|^>=0 (1.10) 
for all physical states. Although Ga{x) do not commute, their commutators 

[Ga{x), G,{y)] = gUAx - y)G,{x) (1.11) 
allow to set them all simultaneously zero. 
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2 Rigid Gauge Rotor. 



In this section I will discuss the rigid gauge "rotations". Classically I define them as 
gauge field configurations of the type 

AUx,t) = U(x,t)aUx)U+(x,t) + -U(x,t)d'U+(x,t) (2.1) 

9 

where a*(x) are t-independent, fixed as far as their x-dependence is concerned, "rigid" 
fields which I do not specify and leave them arbitrary for the moment. Eq. (|2.1|) is the 
simplest example of the transformation to the "body fixed" frame of the local gauge 
symmetry in which I have assumed that the dynamics of the field in this frame is very 
stiff so that the field can be approximately replaced by its static average. In general 
a* is of course dynamical but should be viewed as constrained since U{x, t) already 
contains a third of the degrees of freedom. For non rigid a* there is no obvious choice 
of the body-fixed frame and it can be constrained in a variety of ways, say, = 
(axial gauge), 9ja* = (Coulomb gauge), etc. In our language these different gauge 
fixings correspond to different "rotating" frames. Since they are "non inertial" the 
dynamics will look very differently depending on the choice of the frame and different 
fictitious forces, the analogue of Coriolis and centrifugal forces, will be present. I am 
planning to discuss these issues elsewhere. 

With the anzatz ( p.l| ) the covariant derivatives are 

D' = d'- igA' = U{x, t)d'{x)U+{x, t), (2.2) 

with fixed, rigid 

d\x) = d' -iga\x). (2.3) 

Inserting (|2.1| ) in the Hamiltonian (|1.1|) one finds that the gauge invariant potential 
term Ei,a{^ilx)f ~ j:i,jTr[D\ D^]^ = EijTrld^ d^^ is independent of the U\ 
i.e. it is fixed, nondynamical in this model. The dynamics of the gauge field is 
governed by the kinetic energy, i.e. the term with the electric field in (|1 . 1| ) . Using 
do{UdiU+) = {i/2)U{diu)U+ with u = 2iU+doU one finds 

-E' = doA' = ^U[u;,d']U+, (2.4) 

and therefore the kinetic energy in (|1 . 1|) is 

J I d'xTr{doA'Y = I d'xTr (u;[d\ [d\ to]]) (2.5) 
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where I have disregarded surface terms, ignoring for the moment possible non van- 
ishing fields at infinity, non trivial topologies and other global issues (cf. below). 

The double commutator in ( p^.5[ ) with the summation over all indices, x, i and 
the color is the straightforward generalization of the familiar double vector product 
summed over all particles indices in the moment of inertia tensor appearing in the 
kinetic energy of rigid space rotations of system of particles with fixed relative po- 
sitions. Following this analogy the energy (|2.5|) of the rigid gauge rotations can be 
written 

Erot = (fxTr{ujIuj) (2.6) 
where the moment of inertia is defined as a differential matrix operator such that 



1 



d\[d\uj] 
{d} — ig[dia\ uj] — 2ig[a\ dioo] 



a , [a ,u 



(2.7) 



To obtain the corresponding Hamiltonian one can use the gauge field part of the 
generators (p.. 81). Using Eqs. ( |2.1| ), ( |2.4| ) and the definition ( |2.7| ) one finds 



= [d' - igA\ E' 



29 



-U 



d\ [d\uj]] U+ = -2gU{Iuj)U^ 



(2.8) 



Defining the gauge generators in the rotating frame G = U~^G^U, expressing lj = 
—I~^G/2g from (|2.8| ) and substituting in (|2.6|) one finds the Hamiltonian of the rigid 
gauge rotor 



H, 



A 
rot 



16g 



^ J d^xd^yTrG{x)r\x,y)Giy) = ^J d'xd'yG^ix)I~,\x,y)Gtiy). 



(2.9) 

where I^^{x^y) = (l/4)Tr(A"/ ^(x,?/)A*) is proportional to the inverse of the opera- 
tor —dl^^dli, with dl^jj = diSab — gfabcdl I assumed that this operator does not have 
zero eigenvalues. In a more careful way of handling fields at infinity one should avoid 
the integration by parts in ( p.5| ). The inverse "moment of inertia" operator is then 
replaced by a less transparent 



d^y \dlc{y)ha{y^x) \ \dlAy)h'l'{y^x') 



(2.10) 



Most of the following results remain valid for both forms of this operator. 
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The meaning of the preceeding expressions is quite obvious. They are the field- 
theoretic generahzation of the standard rigid rotor resuhs. The unbroken local gauge 
symmetry of QCD means that there are free SU(N) color gauge "rotations" at every 
space point. Expression (|2.9|) shows that the "rotations" at different points as well as 
around different color axes are coupled via the non diagonal elements of the moment of 
inertia "tensor" Iab{x, y) in the manner similar to the coupling between rigid rotations 
around different space axes in systems of particles. 

It does not seem to be useful to diagonalize the operator in ( |2.9|) . The standard 
diagonal form of the rigid rotor Hamiltonian , i.e., H = ^J^a^l/^a can be usefully 
achieved only in the case of rotations corresponding to a single SU(2) group to which 
the familiar rigid space rotations belong. Diagonalizing the moment of inertia in the 
case of higher groups will introduce combinations of the generators multiplied by ma- 
trices of orthogonal rotations. These in general will not have the group commutation 
relations. Already for a single SU(3) the group 0(8) of rotations in the adjoint space 
needed in order to diagonalize the moment of inertia is much larger than SU(3). 

The actual values of the moment of inertia depend on the rigid configuration a^{x) 
of the gauge field via the expression ( p.7| ). This comprises the "free parameter" of 
the rigid gauge rotor model. For abelian theory or alternatively in the limit g 
the inverse of I{x,y) appearing in Eq.( |2.9| ) is just the Coulomb propagator. In the 
opposite large g or long wavelength limit I{x,y) becomes a local tensor given by 
the last term in ( |2.7] ) which is obviously the SU(N) generalization of the moment of 
inertia expression. 

An important feature of the Hamiltonian ( p.9|) is that despite the severe limitation 



of the allowed gauge field configurations imposed by ( |2.1| ) it remaines gauge invariant. 
This is because ( p.9|) depends on G rather than G^. Under a gauge transformation 
U — > SU, G^ transforms as SG^S^ so that G{x) and therefore H^^^ stay invariant. 
The gauge invariance of ( p.9|) is the simplest illustration of the usefulness of the 
introduction of the body fixed frame. One can freely approximate the dynamics in 
this frame without fears of violating the symmetry with respect to which the frame 
has been defined, i.e. the local gauge symmetry in the present case. 

Consider another transformation, U — >■ US. Referring to Eq. one can in- 
terpret this transformation either as the change of U i.e. the transformation of 
the intrinsic frame with respect to the rigid "shape" a* or as the change of a* , 
a* — >■ Sa^S'^ + ^Sd'^S'^, i.e. the transformation of the intrinsic "shape" with respect 
to the intrinsic frame. Such transformations obviously form a group of local SU(N) 
gauge transformations which I will call the intrinsic or "body fixed" gauge group to 
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distinguish it from the "laboratory" gauge group of the ordinary gauge transforma- 
tions. According to two different interpretations of the intrinsic gauge transformations 
given above one has two options. One option is to regard G's as the generators of the 
intrinsic group. They act on the dynamical variables U but they have a disadvantage 
in that the "laboratory" group is not completely independent of such an intrinsic 
group, e.g. they both have identical Casimir operators. Another option is to formally 
introduce operators which gauge transform the intrinsic variables a*. They will have 
the same form as G^^s but with a* replacing A\ Defined in this way the intrinsic 
group will be completely independent of the "laboratory" gauge group but will act 
on nondynamical variables which do not appear in the wavef unctions. Convenience 
should dictate which one to use. 

The above introduction of the intrinsic vs " laboratory" gauge groups is obviously 
quite general with e.g. the definition of G{x) being independent of the rigid rotor 
restrictions set by fixing a* to be nondynamical in ( ^.11 ). Unlike the local gauge 
symmetry in "laboratory", the symmetry in the "body fixed" frame can be broken. 
E.g., the gauge invariant Hamiltonian (|2.9|) is in general not invariant under the 
transformations of the intrinsic gauge group. This is a simple example of the situation 
to which I referred earlier as a possible existence of "deformation" vs impossibility 
of the symmetry breakdown in the context of non abelian local gauge theory. The 
character of the deformation can be classified using the intrinsic gauge group, e.g. 
in classification of possible "deformed shapes" of the rigid gauge rotor ( |2.9| ) by the 
transformation properties of the moment of inertia Iab{x,y) under this group. Here 
I obviously adopt the second interpretation of the intrinsic group. The invariance 
of Iab{x,y) under all intrinsic transformations would be analogous to the spherical 
rotor limit in the space rotation case. The invariance under a continuous subgroup of 
the intrinsic group is the analog of the axial symmetric rotor, etc. Discrete intrinsic 
subgroups should also be considered. 

Consider a rigid gauge configuration a* which is a gauge transform of a*, a* = 
^(a* + {i/ g)di)S'^ . The Hamiltonian (|2.9| ) will have the same form with the same mo- 
ment of inertia but with G replaced by S~^GS. The eigenvalues of this transformed 
Hamiltonian will not change and will therefore depend only on gauge invariant com- 
binations of the rigid a*, i.e. on the Wilson loop variables TrPexp{ig § a^dxi). 
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3 Static Quarks 



So far I have discussed the rigid gauge rotor hmit of only the first term in ( |1 . 1|) . 
The resulting H^^^ is relevant for the discussion of very heavy quarks. They can be 
considered static as far as their translational motion is concerned. They still have 
a wavefunction describing the motion of their color degrees of freedom. Because of 
( |1.10| ) this motion is coupled to the "rotations" of the gauge field which I will treat 
using ( p^ . 

In the limit of m — oo the quark kinetic energy term q^apq and the quark color 
current coupling q^aX^q in Eq.( |l.l|) can be neglected and the resulting Hamiltonian 
decouples into a part containing the gauge field and another containing the quarks, 
H = Ha + Hg, where is the first term in (|1.1| ) and Hg = m J (Pxq~^{x)Pq{x). The 
coupling appears only via the Gauss law constraint, Eq.( p..lO| ). The wave function 
can not be taken as a product \I' = \I'(g)\E'(y4) but should be a local color singlet. In 
the representation in which 




qa{x) = aa{x) ^ 
with 

{aa{x),a'^{y)} = 5o,p5{x - y) 

{ha{x),h']^{y)} = 5cp5{x -y),etc... 

the eigenfunctions of Hq are trivially written down. Consider, e.g., 

\m{q) >= \vac{q) > = |0 >, (3.1) 
|^(g) >= |xo,a > = a+(xo)|0> (3.2) 
\^{q)>=\xo,a;yo,(3> = a+(xo)6^(yo)|0 > (3.3) 

These wave functions describe respectively zero quarks, one static quark at Xq with 
color component a and a static quark - antiquark pair at Xq and yQ. It is easy to form 
local color singlets with these wave functions. For e.g. the quark- antiquark pair it is 

1^' >=J2<i/,,,a;yoAA)\xo,a;yo,P> (3.4) 

0/3 



11 



with the wave functional '^xo,a;yo,i3{A) satisfying 

Gtix)^x,,^.,y,AA) = g5{x - xo)^v^,„,„^,„,^(A) + g5{x - yo)^vl/^^^^^^^_^, (AX3.5) 

where A^^ = — A^^. The wave functional of the gauge field should be a singlet at 
every point in space except at the position of the quarks where it should transform as 
N and N multiplets of SU(N). This constraint together with the Schrodinger equation 

HA'^.o^o.^yoAA) = E^,,,o.;yoAA) (3-6) 

completely defines the problem for the gauge field. 

In the rigid gauge rotor limit Hj^ is given by Eq. (|2.9|) . The wavefunctions of 
this Hamiltonian are general functionals \E'[?7(x)] of the SU(N) matrices Ua/six). 
Their scalar product is determined by functional integration over the t/'s with the 
corresponding group invariant measure. The vacuum wave functional must obey 
G^{x)'^^ac[U] = 0. This means that it is a constant independent of Ua/six). Since 
also Ga{x)'^vac[U] = the vacuum energy is zero according to (p.9|) . Regarding the 
parametrization of the t/'s in terms of the appropriate Euler angles of the SU(N) rota- 
tions at every space point, the constant \Ef„ac[f^(^)] nieans that all the "orientations" 
of U (x) at all points are equally probable, i.e. there are no correlations between the 
"orientations" of the rigid gauge rotor at different points. This is as "random" as the 
distribution of the f/'s can get. The absence of correlations is the property only of 
the vacuum. For other states the "orientations" of the gauge fields at different space 
points are correlated via the "moment of inertia" operator. 

In order to discuss the wave functions with non zero number of quarks it is suffi- 
cient to know some simple properties of the gauge generators G^{x) and G'a(x). Since 
Ga{x) are gauge scalars, they commute with G^{x), 

[G^{x),G,{y)] = (3.7) 

which means that together the generators of the " laboratory" and the intrinsic gauge 
groups provide a complete set of commuting quantum numbers for the wave func- 
tionals "^[Ua/six)]. Indeed, since the Casimir operators for G'^'s and G's coincide one 
has e.g., for the SU(2) the {G'^{x)Y, G^{x) and Gsi^x), i.e. three local commuting 
operator fields for the three fields of the Euler angles needed to specify the U{x). 
In the SU(3) one has eight fields of the "Euler angles" and eight local commuting 
generators made off G^'s and G's - the two group Casimir operators, one Casimir 
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operator of an SU(2) subgroup for G s, say Z]a=i(^a )^(^) ^^e corresponding 
one for the G's and respectively two pairs of the Cartan generators - Gf{x), Gf{x) 
and 6*3 (x) and Gs{x). This counting continues correctly for any N, i.e. — 1 for the 
SU(N) Casimir operators, 2((A^ — 1) + (A^ — 2) + ... + 1) for the Casimir operators 
of pairs of SU(A^ — 1)...SU(2) subgroups and 2(A^ — 1) for the Cartan generators. 
Altogether there are A^^ — 1 local commuting operators as needed. An eigenfunction 
of this complete set of operators is the Wigner function D^j^,{U{x)) of U at a cer- 
tain space point. K and K' are the quantum numbers of the "laboratory" and the 
intrinsic groups and L determines the representation. 

Since under an infinitesimal gauge transformations f/ — (1 + iea{x)^)U one can 
easily verify that 



[Ga{x),Uo.M] 

Ga{x),U^p{y) 

Ga{x), Uap{y) 

Gaix),U^^iy) 



A" 

g^ix - y)^U^/3{x) 



-96ix-y)U^^ix)^ 



X"' 

g6{x - y)Uayix)-^ 
A" 

-g6{x ~y)^U+p{x) 



(3.8) 



All the operators in the rigid gauge rotor model are functions of the G"s and ?7's. 
E.g. consider the electric field operator. According to Eg. (^.41) it is 
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U[d\I-^G]U^ 



(3.9) 



where I have expressed u in terms of G using (|2^ 

Using the relations ( |3.8| ) it is easy to write the general form of the wave functions 
for a single quark and for a quark-ant iquark pair, 



f^[U] = Uay{xo)U^f3{yo)Cj5 



(3.10) 



They satisfy the conditions ( p.5|) following from the Gauss law with constant coef- 
ficients and c^s which give the probability amplitudes of the intrinsic quantum 
numbers 7 and 6. They should be normalized, I] |c^P = 1; ^75^ = 1 to assure the 
normalization J d[U{x)]\'^[U]\'^ = 1 
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These amplitudes must be found by solving the corresponding Schrodinger equa- 
tions but before describing this I wish to remark that the above form of the wave 
functions is valid also when the limitation of the rigid gauge rotations is relaxed and 



the most general gauge configurations are allowed. The parametrization (|2.1| ) is still 
very useful but now with fully dynamical fields a* the variation of which should be 
limited only by a "gauge fixing" condition to avoid overcounting as described above. 
The dynamics will of course be that of the full QCD but the wave functions of the 
static quark and the quark-antiquark pair will have the same form ( p.lO|) . The dif- 
ference will be that amplitudes and c^s will be functionals of a*(x) describing 
the space and color fluctuations of the "string" attached to the quark or between 
the quark and the antiquark. In the rigid gauge rotation case there are only color 
fluctuations described by constant amplitudes. 

For quarks in higher representations the wave functions have the same form with 
U replaced by the appropriate Wigner D-function. E.g. in the adjoint representation 

^xo,a[f/] = Tr(f/(a;o)A'^f/+(xo)A')c6, 

etc. 

I will now derive the Schrodinger equation for the string amplitudes and c^s. 
Acting with the Hamiltonian (2^) on ( 3.101) , using (|3.8|) and the orthogonality of t/'s 



with respect to the integration over the group, / dUU*f^U^u = Sa^j.5i3u, I find 

Qa-y (3^0) '-'7 E^ai 

where I denoted 

Q„^(a;o) = ^/-,i(a;o,xo)(A"A^),^, (3.12) 

-P«/3,7/i(a;o,l/o) = 2-^ab^(^o,Z/o)A[5,^A|^^. (3.13) 

In SU(2) Qa-y takes a particularly simple diagonal form, Qa-y = 6a-y{l/4:)I~J^{xo,xo). 
and is the eigenvalue for a single quark. For quarks in e.g. adjoint representation the 
lambda matrices in the expressions above are replaced by the corresponding group 
generators ifabc- The first two terms in the second line of ( p. 11 ) are the quark and 



the antiquark self energies whereas the last term is their interaction. In QCD one 
expects that terms like Q are inflicted by the long and short distance divergences and 
should be properly regularized which I will assume for the rest of the paper. I will 
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further assume the translational invariance of Q, i.e. its independence of xq. One 
can then rewrite Eg. ( p. Ill) by transforming it to the basis in which Q is diagonal. 
Defining its eigenvectors Qfo^") = enfo*^"^ and expanding c^p = dmub'^^^b^^^* one finds 

< kl\P{xo,yo)\mn > dmn = {E — — e^)dki (no sum over k and 1) (3.14) 



where < kl\P\mn >= P„^,T,^6(™)65r)*6(f)6{j^*. The Schrodinger equation ( p.l4| ) is 



jV^ X A^^ matrix equation and the most interesting question of course concerns the 
dependence of its eigenvalues on the distance |xo — yo\ for various possible choices 
of the rigid gauge field configuration a*(x) on which the matrix P depends. I will 
address this question in the next section. 



4 Choices of The Rigid Field. Mean Field Equa- 
tions. 

The rigid configuration if it exists in QCD must reflect the properties of the gluon 
condensate of the vacuum. One of the more accepted views of the QCD vacuum is 
that this is a condensate of non trivial topological configurations - the Z(N) vortices, 



c.f.,|]ll|. Although such configurations are easily incorporated in the above formalism 
I was not able to overcome technical difficulties in working out a theory of their 
condensation. 

On a heuristic level each Z(N) vortex carries a unit of flux of the colormagnetic 
field. Condensation of the vortices presumably means that there is a non zero average 
of this field in the vacuum. Of course due to unbroken local gauge symmetry it must 
undergo free "gauge rotations" at each space point. In the ground state this means 
that there are equal probabilities of all the "orientations" yielding zero average value 
in the laboratory. The finite average value of the condensate field can only be " seen" 
in the "body fixed" frame and should appear in this picture in the manner similar to 
a* in the expression ( [2.1D for our rigid gauge rotations. The field strength 

B\x) = U{x)b\x)U+{x), U = ^-e,jk[d\ d% (4.1) 

also averages to zero in the ground state but has a non zero value 6* in the "body 
fixed" frame. 

Via the dynamics of U{x) the anzatz (p.l|) leads to colorelectric field (|3.9| ) which 



propagates away from points where G{x) is non zero, i.e. from the location of static 
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quarks. The propagator of this field is controlled by the condensate field a* which 
enters the expressions for and d^. This propagator is a long range Coulomb 
potential for zero a* and is a Gaussian for a* corresponding to a uniform colormagnetic 
6* (cf., below). The screening of the propagation range of the colorelectric field 
in the presence of the colormagnetic "condensate" 6* is reminiscent of the dual to 
the Meissner effect of screening of a magnetic field by the electric condensate of a 
superconductor. This possibility of the dual Meissner effect is of course a standard 
scenario for confinement in QCD. It is expected that tubes of fiux of the colorelectric 
field are formed which connect quarks and make their energy depend linearly on the 
distance. 

In the present formalism a way to attempt to model the formation of a confining 
string is to look for such a configuration of the rigid field a*(x) for which the prop- 
agator behaves roughly speaking as one dimensional for large separations along 
some given line in space at the end of which quarks can be placed. This means that 
a sort of magnetic "wave guide" should be constructed so that the Green's function 
of the operator -d^^ = -dl^d^f, = -{diSac - gfacc'al,){di - gfcbvay) is asymptotically 
oc |a; — x'l along, say, one of the coordinate axes. In order to see the difficulties in 
finding such a configuration consider for simplicity 2 space dimensions and choose 
a} = c{y) and = with an arbitrary c{y). This choice corresponds to the color- 
magnetic field b{y) = dyc{y) depending only on one coordinate y. The operator to 
invert is then 

- (d. - ^gciy) . - (4.2) 

where I denoted the color spin matrices Fj^^ = ifbac and c{y) ■ F = Ca{y)F"-. The 
propagator is then 

where Xn{k,x) and e„(fc) are solutions of 

l-d'y + {k + c{y) ■ Ff]xn{k, y) = e^{k)xn{k, y). (4.4) 

In order to achieve the desired confining behaviour of the propagator the sum in (^4.3|) 
must be ~ k'^ for — > 0. The simplest is to assume that the lowest eigenvalue of 
([4.4| ), eo(A;) should vanish as fc^ for small k. However the operator in (f4.4|) is a sum of 
squares and does not have zero eigenvalues for non trivial regular c{y). It is also not 
symmetric in k for small values of k but this seems to be less of a problem. The same 
conclusions seem to hold in 3 space dimensions. It is quite possible that perhaps a 
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singular configuration a' exists which leads to zero eigenvalue in ( ^.41 ) at zero k but 
I was not able to find it. 

The strong coupling limit of lattice QCD suggests that quarks in the fundamental 
reperesentation are confined whereas they are only screened if put in the adjoint 
representation. This crucial difference comes from fairly simple quantum mechanics 
of color degrees of freedom related to matching of group representations in neighboring 
lattice points. In our rigid gauge rotor model a similar simple quantum mechanics of 
colors is retained. As a result the eigenenergies of a systems of static quarks will be 
determined by different combinations of the color components of depending on the 
representation of the quarks. E.g., as already mentioned in Section 3 when the quarks 
are taken in the adjoint representation the A matrices in the expressions for P and Q 
in the Schrodinger equations ( |3.11| ) are replaced by their adjoint counterpartners F. 

In order to find the optimal a* in a systematic way one can follow a variational 
approach and minimize the ground state energy of the rigid gauge rotations for fixed 
positions of static quarks. This energy is given by a sum of the lowest eigenenergy of 
H^g^, Eq.( |2.9|) and the colormagnetic energy given by the second term in ( |1.1|) with 
rigidly "rotating" A{x), Eq.dg^), i.e. 

1 

Variation of this expression gives 



E[a'] = Erot[a'] - / d?xTr[d\ d^f (4.5) 



dj^^-ig[a\n = \^, (4.6) 

where f^^ = {i/g)[d\d^] . Eq.(^]B|) is obviously gauge invariant. In the vacuum Erot 
is zero and the minimization of the second term simply gives the classical equation 
for a* in the vacuum. For a quark- antiquark system Erot is non trivial and depends 
on the distance between the quarks. I plan to discuss the solutions of the equation 
( [4. 61 ) and their relation to confinement elsewhere. 

In the rest of this Section as an illustration of a simple choice for the rigid field 
a* which allows to obtain some analytic results I consider it to be diagonal, al^p{x) = 
SajjCiai^)- The moment of inertia operator with such a* is 



1 



1 



Using Green's function satisfying 

{d' - ig{d^ - a'p)f J„^(x, y) = -5{x - y), (4.8) 
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and following the procedure leading to Eq. 
Hamiltonian in this case 



one finds the rigid gauge rotor 



H = ^ I d^x<fyGo,f3{x)Ja(3{x,y)Gpc,{y)- 



(4.9) 



The Schrodinger equation for the static quark-antiquark wave function ( p.lOl ) has the 
form ( p. 14 ) with 



Jal3{XQ, Xq) — — \ 2Jaa{XQ, ^o) — J^'^ JlS/siXo, Xq) 



(4.10) 



{xo,yo) 



(4.11) 



1 



The diagonal components of J are simple Coulomb propagators independent of the 
color so that the expressions for Q and P can be simplified further but I will not go 
into the details of this. Instead I will now consider the choice of a* which corresponds 
to a much discussed in the literature situation of a uniform colormagnetic field. I 
emphasize that in the present model this field is uniform in the intrinsic, "body 
fixed" frame. For simplicity I will first work in 2+1 dimensions and will try to extend 
to 3+1 in the next section. I set 

1 



2 ^« ^ijX'^ 



(4.12) 



where the space indices i,j presently run over the values 1 and 2. In two space 
dimensions one can take h diagonal in color since the transformation diagonalizing it 
is a part of ?7's in ( p.lj) . Explicit expression for J is easily obtained in this case from 
the known Green's function of a Schrodinger equation in a constant magnetic field, 
cf. Ref.0, 



Jap{x,y) 



1 

— f 

47r 



°° ds 



sinh s 



-(l9ba/3l/4){a;-j/)2 coths 



(4.13) 



where b. 



'af3 



hp. For X = y this expression is independent of color indices. It 



must be regulated to prevent the divergence, e.g. 

g'^N ds 



27T 



sinh s ' 



(4.14) 
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where sq is a regularization cutoff. Although J(x, y) does not depend only on 
the distance \x — y\ the Schrodinger equation (|3.14|) with P and Q based on such 
J is translationally invariant. Shifting the coordinates by say a vector h and 

simultaneously performing a gauge transformation of the wave function c^^ — ^ 

J 



exp 



i{gba/'2)eijh\xQ — yl) leaves Eq. (|3.14| ) invariant. The integral in the expres- 



sion for J{x,y) can be expressed in terms of the Bessel function Ko{\g{ba — hp\z^ /A) 
with z = X — y and for \g{ha — bi3)\z^ — oo it has the following asymptotic form 



1 



4v^ 



{2\g{h^ - 6/3)1)-'/' exp \-\g{h^ - hp)\z^/A 



(4.15) 



For finite values oi g\ba — hp\ it decreases as a Gaussian at large separations z. This 
should lead to a similar decrease of the eigenvalues of ( 3.14|) - an entirely unsatisfac- 
tory behavior as far as the confinement is concerned. In the next Section it will be 
seen that the situation may be different in the large N limit. 



5 Large N Random Color magnetic Fields. 

As mentioned in the Introduction rigid gauge field rotations should be relevant for 
QCD in the large N limit where it is expected that a master field configuration 
dominates the vacuum, |Q. As in the case of a condensate such a configuration can 
not be just some fixed gauge field potential ^^(a;). It must be allowed to undergo free 
gauge "rotations" exactly as a' in Eq. (|2.1|) since the gauge invariance is not expected 
to be broken in the large N limit. The dynamics of these rotations can not be "frozen" 
and must be described by the gauge rotor Hamiltonian considered in Section 2. These 
"rotations" induce an interaction between quarks as was shown in Section 3 for static 
quarks and will be demonstrated for dynamical quarks in Section 6 below where it 
will also be shown that in addition a* appears as a background field in the Dirac 
operator. 

Another important consideration is that for large N there is a large number of 
degrees of freedom operating at each space point which introduces statistical elements 



in the theory, cf. Refs. ||T9|, |2^, Experience with this limit for simple systems 
indicates that two types of gauge invariant physical operators should exist, analogoes 
to macroscopic and microscopic observables in thermodynamics. The former depend 
on finite (relative to N) number of dynamical variables and involve sums over all 
labels of the degrees of freedom, i.e. the color indices. A simple example is aj^^a^^, 
etc. Operators without such summations, e.g., a^^ with fixed a and /? belong to the 
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second type which must be regarded as microscopic observables hke ,e.g., a coordinate 
of a particle or a single spin variable in thermodynamic systems. The fluctuations of 
the macroscopic operators are suppressed and expectations of their products factorize 
at = oo. This is not so for microscopic observables. 

On the basis of these considerations one can adopt the following point of view. 
After allowing for free gauge rotations according to (|2.1| ), i.e. after transformation 
to the body- fixed frame, one should consider a^{x) as static random matrix functions 
described by a probability distribution P[a\x)]. This distribution can be determined 
following the ideas of the random matrix theory, cf . Ref . |^ . To this end one should 
introduce the amount of information (negative entropy) 

/{P[a]} = J Dfi[a\x)]P[a\x)]lnP[a\x)] (5.1) 

associated with the P[a^{x)]. Minimizing /{P[a]} subject to suitably chosen con- 
straints on macroscopic-like variables should determine the least biased distribution 
P[a*(x)]. As in statistical mechanics the large N factorization should then simply 
appear as a consequence of the central limit theorem. 

There are two crucial questions which need to be answered in following this proce- 
dure - what is the appropriate measure in the integral (|5.1| ) and what are the variables 
which should be constrained. I hope to address the general answer to these questions 
in the future work. Presently I will illustrate how the procedure can be put to work 
for a uniform colormagnetic field, Eq. ( |4.12| ). 

In the limit of large N only the first terms in the expressions for P and Q above 
should be retained and the Schrodinger equation ( p.l4| ) for diagonal components of 
the string amplitude becomes 

- 2(7^ ^ j^^(xo, yo)cpp = {E- El)c^^, (5.2) 

where E^ = Ig^Y.^Jani.XQ.Xo). The non diagonal string amplitudes decouple and 
satisfy a trivial equation (i?° -|- E^)Cfj,i, = 2Ec^u the eigenvalues of which are simply 
the sums of the selfenergies. Without careful treatment of long and short distance 
regularization in the large N limit one can not reliably discuss these eigenvalues 
taken separately and I will concentrate on Eq.( ^.2|) . Using translational invariance 
and writing this equation for xq = and yo = z one obtains 

E/ r^^"^"^'""''^'^'^^'^°''^^/^/^ = -TF2^- (5-3) 
rrl "'0 47r smh s 2q^ 



20 



as the large N limit of the Schrodinger equation for a static quark - antiquark pair in 
the rigid gauge configuration corresponding to a uniform colormagnetic field in 2+1 
dimensions. One must still specify the large N scaling of various quantities which 
enter this equation. Provided each term in the sum on the left hand side is of the 
same order of magnitude I get the standard scaling of the coupling constant requiring 
that g'^N is held fixed. In the exponential of the integrand one can then extract the 
finite combination g = gy/N. The problem is then to determine the scaling and 
in general the entire distribution of the field components ba/VN. Regarding the 
behavior at large separations z one notes that if the limit of ^ cxd is taken first in 
such a way that the differences \ba — hp\/\/N decrease then the Gaussian decay can 
possibly be prevented. 

I use this example to demonstrate how the ideas about the statistical nature of the 
large N limit can be used to determine the distribution of the components ba/y/N- 
I consider what happens with the Wilson loop W{C) = jj < TrP exp{ig §q A^dxi) > 
in the present theory. Choosing the loop perpendicular to the time axis, inserting 
( |2.1|) in W{C), using its gauge invariance and the explicit form (|4.12| ) of a* one finds 



1 rao 

W{C) = —Y,e'3^-^= dbp{b)e''^^ (5.4) 

1\ ^ J —oo 

where 5* is the area of the loop and p{b) = {l/N)Y,a=i^ib — ba) is the density of 
the field components. In the large N limit p{b) can be approximated by a smooth 
function provided the range of variations of b does not grow with N. Assuming this 
one easily finds simple expressions for p{b) , e.g. Lorenzian 



which lead to the area law dependence of the Wilson loop, W{C) = exp{—gboS). The 
combination gb^ plays the role of the string constant. The placing of A^'s in ( |5.5|) was 
chosen in such a way as to have this constant finite for A^ ^ oo. 

The choice ( |5.5|) is the simplest possible. Any meromorfic function p{b) with poles 
in the upper plane will give the area law with the string tension controlled by the 
position of the pole closest to the real axis. One can also take functions with other 
type of singularities in the upper complex plane, etc. The simple choice ( ^.5|) gives 
area law for any S, missing entirely the asymptotic freedom behavior at small S. One 
can attempt to correct this by choosing more involved expressions for p. A much more 
serious problem is that the space oriented Wilson loop may not be a good measure 
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of confining properties in the model where one has never worried about the Lorenz 
invariance. 

Adopting any form of the "single component" density p(6) still leaves the distri- 
bution of the values of ba needed in, e.g., Eq. (|5.3| ) largely undetermined. Using 
statistical concepts described in the beginning of this Section one should view b^s as 
random quantities and introduce their joint probability distribution P{bi,b2, ■■■,b]sr) 
which should be such that p{b) is reproduced but is otherwise maximally random, i.e. 
contains least amount of information. The question immediately arises as to whether 
p{b) is the only quantity which should constrain P{bi, ...,bN) and what is the com- 
plete set of such constraints. In the absence of general answers I take p{b) controlling 
W{C) as an example and determine the distribution P{bi, bj^i) by minimizing the 
appropriate negative entropy (information) with this constraint. 

The quantities b^s are eigenvalues of a hermitian, in general complex matrix. The 
information content of a probability distribution P(&i, b]\f) of such eigenvalues is a 
well studied question, cf. Ref . pl| . It is 



rf/i[6]P(6i,...,M^^m,-,M (5.6) 

where the measure is dp[b] = const na>/3 — bf3\'^dbidb2...dbN, reflecting the re- 
pulsion of the eigenvalues. Minimizing ( ^.6|) under the condition of a given p{b) = 
{l/N)Y.aip-b^) one finds 

( N N ^ \ 

P{bi,...,bN) = canst exp\Y,Mba-bp\-2NY, ln\ba ~ b'\p{b')db' \ . (5.7) 

Using, e.g. Eq. ( |5.5| ) for p(6) this expression becomes explicitly 

(N N \ 

Y,ln\ba-bf3\- NY,Hbl + Nbl)\ . (5.8) 

The constant in front of this expression must assure the normalization of P and can 
be calculated by the methods described in Ref.p^]. Using the standard interpretation 



of P{bi, ...,6Ar) as a partition function of a fictitious Coulomb gas one can say that 
the "particles" ba are "repelled" from each other by the first term in its exponential 
but are kept within the interval bo by the second term representing the interaction 
with the background "charge" distributed according to p{b). The average distance 
\ba — bpl ^ bo/N becomes very small in the large N limit. The Schrodinger equation 



22 



( |5.3| ) is now a random matrix equation with the probabihty distribution of its elements 
controlled by the P{bi, Bn) above. The actual numerical solution of this equation 
is now in progress. 

In a similar way one can consider rigid gauge configuration representing uniform 
colormagnetic field in 3 + 1 dimensions. This field in the intrinsic frame corresponds 
to the choice 

a'{x) = -eijkVx^ (5.9) 

where now however the three color matrices U in general cannot be assumed diagonal. 
For such non-diagonal colormagnetic field the inversion of the moment of inertia 
operator —{l/Ag'^) [d*, [rf*, uj\\ requires a solution of a matrix differential equation. This 
equation simplifies considerably if 6*'s are nonetheless restricted to be diagonal, = 
^a/sba- Then the equation ( |4.7] ) is still valid with the index i now running from 1 to 
3. In the following equation (|4.8| ) for the Green's function one should just replace 
{ba ~ bp)eijX^ by eijk{bi^ — h'p)x^ . The expression for this Green's function is known 
and one can repeat all the steps leading to the static quark-antiquark Schrodinger 
equation which is the analog of Eq.( p.3| ) in 3+1 dimensions. 

Turning again to the Wilson loop one finds in this case f^a^dxi = VSj with 
Sj = {l/2)ejki §cx'^dxi so that {§^a^dxi)^ = baScosOa (no sum over a) where S = 



J2i{S''y, ba = yHiib'^aY and 9a - the angle between the vectors and S"" at a given 
a. S is the area of the loop when it is planar and is related to the minimal area in 
general. The Wilson loop is 

W{C) = — < g'ff'S'fea cos da ^_ 



j^Y.'^^ f d{cos9a)e''^^''"'°''^" = (5.10) 

~^ I ir/'Wsin(^&^), 
gb Jo b 



where p{b) = (l/N) J2a^{b — ba) is the density of the positive lengths of the color 
components of the vector b\ In (|5.10| ) I have performed the angle averaging which 
must be present in the vacuum wavefunction. One can easily choose p{b), e.g. the 
square of Lorentzian 

^^^^ = .{b^ + Nbiy ^'-''^ 
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which gives the area law W{C) = An exp{—gboS). This choice is again not unique and 
gives the area law for any 5". It has a powerlike tale as opposed to the perturbative 
Gaussian. 

The statistical arguments for finding the entire distribution of 6^ can be used in 
3 + 1 dimensional well with the difference that in this case the density of the 

lengths of the vectors b]^ is fixed by, e.g. Eq. (|5.11| ) and their directions are distributed 
isotropically. 



6 Dynamic Quarks 

Dynamic quarks can be easily included in the rigid gauge rotation model. For this 
I define quark fields in the "rotating frame", q = Uq, use Eq. in the second 

term of the QCD Hamiltonian ( |1 . 1|) and replace the first term in it by (|2.9| ). Using 
moreover the Gauss law constraint ( |1.1CI| ) I can write the original QCD Hamiltonian 



1.1) in the rigid gauge rotor limit as expressed in terms of the quark fields only. 



Hrot = ^Jd^xd^ypa{x)lJ{x,y, [a'])pb{y) + J d^xq+{x)[a'{p' - ga\x)) + (3m]q{x), 

(6.1) 

where pa = q^X^q are the color quark densities in the rotating frame. The Hamil- 
tonian Hrot describes quarks with gauge strings attached to them, i.e. q{x) are 
multiplied by U^{x). They move in an external colormagnetic field described by the 
vector potential a*(x) and interact via an instantaneous interaction I~h{x, y, [a*]) also 
depending on a* via Eq. ( p.7| ). The simultaneous appearance of the rigid gauge field 
configuration both as a background and as "inducing" the quark-quark interaction is 
ultimately a consequence of the gauge invariance which requires that non dynamical 



rigid gauge fields appear only in the form (|2.1|) 



Hrot is gauge invariant since the operators q are. Moreover this Hamiltonian 
should only be used in the color singlet sector of the theory since I have used the 
Gauss law to derive it. 

For the vanishing a* the Hamiltonian Hrot describes free quarks interacting 
Coulombically. Also for a general non zero a* Hrot should be regarded as the QCD ana- 
logue of the QED Hamiltonian in which only the instantaneous Coulomb interaction 
between the charges has been retained. Indeed the analogue of the rigid gauge "ro- 
tations" ( pTTD in QED is A^{x,t) = a*(x) ^) with abelian U = exp{igx{x^t)), 
fixed rigid a*(x) and dynamical t). Repeating the steps leading to ( |6.1D one will 
derive in the QED case the Coulomb interaction between the charge densities. 
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Regarding the possible role of a*(x) as the master field in the large N limit one 
has in Hrot a way in which this field should enter the quark sector of the theory, i.e. 
serving both as a background field and perhaps somewhat surprisingly also controlling 
the quark interaction. Following the developments of Section 5 this field should be 
regarded as random. The appearance of a random interaction between the quarks 
means that a possible mechanism for confinement of dynamical quarks in the large 
N limit could be related to the localization of their relative distances. The possible 
connection between confinement and localization has already been mentioned in the 
past but usually in the context of a random background field and not with random 
interactions as appear in the present model. 

The Hamiltonian (|6.1| ) takes exact account of the gauge symmetry. One must 
however also worry about global symmetries. For any N the Hamiltonian Hrot may 
serve as a possible basis for various phenomenological developments. Both for this 
matter and conceptually one must face the issue that allowing for an arbitrary x- 
dependence of various color components of a* in ( |6.1| ) leads to breaking of important 
symmetries such as translational, rotational, Lorenz, time reversal, and various dis- 
crete symmetries. Of course the breaking of continuous space symmetries is not 
uncommon in phenomenology, e.g. the bag model, the quark potential model, the 
Skyrme model, etc. Symmetries can be restored by considering all configurations 
translated by the symmetry and integrating over them using collective coordinates. 
This of course applies to both continuous and discrete symmetries. In the absence of 
the guidance from the symmetries a more dynamical criterion for fixing a* seems to 
be the condition of lowest energy. This leads naturally to a generalization of the vari- 
ational approach of Section 4 in which the variational energy should be replaced by 
the ground state energy EQ[a^{x)] of (the suitably regularized) Hj-ot found for a given 
a*(x). Should the solution a* break a global symmetry, the symmetry "images" of 
this a* will also be solutions and one should "sum" over all of them in a standard way 
thereby restoring the symmetry. This variational approach may be combined with 
the Hartree-Fock method which should allow to calculate EQ[a^{x)\ approximately. 
The Hartree-Fock approximation for fermions was shown to be consistent with the 



large N approximation, c.f., Ref . |P2|| . In a combined approach one should form for 
fixed a* an expectation value of Hrot with respect to a trial state of a chosen color 
singlet configuration of quarks (e.g., vacuum, baryon, etc) which must be a product 
state, i.e. such that the expectation values with respect to this state have a non 
interacting factorized form,i.e., 

< Pa(x)pb(i/) > = KA{< g+(x)g/3(x) >< q^{y)qs{y) > - 
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- < QUx)Qs{y) >< q^{y)qfi{x) >). (6.2) 

For a global color singlet state < >= 5app{y,x), with p{x,y) = -^J^-y < 

q^{y)q-y{x) > and therefore 

< paix)pb{y) >= -26abpix,y)p{y,x) (6.3) 

Following the standard Hartree-Fock routine, cf., Ref.|3^, the single quark den- 
sity matrix can be expanded in terms of a complete set of functions, i.e., p{x,y) = 
J2n fn4'n{x)ipn{y) ■ At this stagc it is customary to add the so called Slater determi- 
nant condition which means that the single quark states ipn have sharp occupations 
/„ = or 1. This condition and its compatibility with the large N limit and Lorenz 
invariance were discussed in Ref. [^. Adopting this condition, using ( |6.3| ) in forming 
the expectation Ehf[p{x, y), a*(x)] =< Hrot > and varying with the respect to ?/'„(x)'s 
with constraints on their normalization one obtains the Hartree-Fock equation 

[ay + (3m]iPnix) -2 J d^yl^a^x, y, [a^])p(x, y)My) = e„V^„(x), (6.4) 

which appears here as a selfconsistent Dirac equation for the quark wave functions 
ipn{x)- Note that in this equation a' disappeared from the Dirac operator and enters 
only the interaction I~^{x,y). Solutions of (|6.4|) determine the optimum p and thus 
Ehf for a given a*. 

Hartree-Fock equation similar to ( |6.4|) has been investigated in the 1-1-1 dimen- 
sional QCD, |]2T]] although there for obvious reasons the field a* was absent and 
I~^{x,y) was simply \x — y\. Both the t'Hooft meson spectrum and the baryon soli- 
ton solutions were found in this approach in 1-1-1 dimensions. For small quark masses 
the baryon was the realization of the skyrmion described in the quark language. If 
successful the approach based on Eq.( |6.4|) can possibly lead to similar results in 2-1-1 
and 3-1-1 dimensions. 

Acknowledgement 

Useful discussions with S. Elizur, S. Finkelstein, J. Goldstone, K. Johnson, A. Ker- 
man, F. Lenz, M. Milgrom, J. Negele, J. Polonyi, E. Rabinovici, S. Solomon and H. 
Weidenmiiller are acknowledged gratefully. Special acknowledgement is to A. Schwim- 
mer for patient teaching and answering my questions and for critical reading of the 
manuscript. Parts of this work were done during visits in Max-Planck-Institut fur 
Kernphysik at Heidelberg and I wish to thank H. Weidenmiiller for warm hospitality 
there. 



26 



References 



[1] A.A. Belavin, A.M. Polyakov, A.Z. Schwartz and Yu.S. Tyupkin, Phys. Lett. 

B59 (1975) 85; G. t'Hooft, Phys. Rev. Lett. 37 (1976) 8; Phys. Rev. D14 (1976) 
3432; A.M. Polyakov, Nucl. Phys. B120 (1977) 429; C. Callan, R. Dashen and 

D. Gross, Phys. Lett. 63B (1976) 334; 66B (1977) 375; Phys. Rev. D17 (1978) 
2717; D19 (1979) 1826; S. Coleman, The Uses of Instantons, Erice Lectures, 
1977; E.V. Shuryak, Nucl. Phys. B203 (1982) 93,116,140; B214 (1983) 237; D. 
I. Diakonov and V. Yu. Pctrov, Nucl. Phys. B245 (1984) 259, B272 (1986) 457, 
B306 (1988) 809, B323 (1989) 53. 

t'Hooft, Nucl. Phys. B72 (1974) 451; B75 (1974) 461. 

E. Brezin, C. Itzykson, G. Parisi and J.B. Zuber, Commun. Math. Phys. 59 
(1978) 35. 

E. Witten, Nucl. Phys. B160 (1979) 57; S. Coleman, 1/N, Erice Lectures, 1979, 
in S. Coleman, Aspects of Symmetry, Cambridge University Press, 1985. 

Yu.M. Makccnko and A.A. Migdal, Phys. Lett. 88B (1979) 135; 89B (1980) 
437(E); 97B (1980) 253; Nucl. Phys. B188 (1981) 268. 

T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 106; G. Bhanot, U. Geller 
and H. Neuberger, Phys. Lett. 113B (1982) 47; 115B (1982) 237; V.A.Kazakov 
and A.A. Migdal Phys. Lett. 116B (1982) 423; A.A. Migdal, Phys. Lett. 116B 
(1982) 425. 

K. Wilson, Phys. Rev. D14 (1974) 2455; in Methods in Field Theory ed. R. 
Balian and J. Zinn-Justin, 1981, p. 261; M. Creutz, Phys. Rev. D21 (1980) 
2308. 

V.A. Kazakov and A.A. Migdal, preprint LPTENS-92/15 (1992). 

J. Kogut and L. Susskind, Phys. Rev. Dll (1975) 395; J. Kogut, Rev. Mod. 
Phys. 51 (1979) 659. 

A.M.Polyakov, Phys. Lett. 59B (1975) 82; Nucl. Phys. B120 (1977) 429; A.M. 
Polyakov, Gauge Fields and Strings, Harwood, New York, 1987. 

G. t'Hooft, Nucl. Phys. B138 (1978) 1; B153 (1979) 141; B238 (1981) 406. 



27 



[12] M.A. Shifman, A.I. Vainstein and V.I. Zakharov, Nucl. Phys. B147 (1979) 385, 
448, 519. 

[13] H.B. Nielsen, Phys. Lett. 80B (1978) 133; H.B.Nielsen and P.Olesen, Nucl. 
Phys. B144 (1978) 376; B160 (1979) 380; P. Olesen, Phys. Scripta, 23 (1981) 
1000, Nucl. Phys. B200[FS4] (1982) 381. 

[14] H.C. Pauh and S.J. Brodsky, Phys. Rev. D32 (1985) 1993; D32 (1985) 2001; 
K. Hornbostel, S.J. Brodsky and H.C. Pauh, Phys. Rev. D41 (1990) 3814; 
R.J. Perry, A. Harindranat and K. Wilson, Phys. Rev. Lett. 65 (1990) 2959; 
S. Brodsky, H.C. Pauh, G. McCarter and S. Pinsky, "The Challenge of Light- 
Cone Quantization of Gauge Field Theory" , OHSTPY-HEP-T-92-005; F. Lenz, 
S. Levit, M. Thies and K. Yazaki, Ann. Phys. 208 (1991) 1. 

[15] K.Johnson, L. Lellouch and J. Polonyi, Nucl. Phys. B367 (1991) 675. 

[16] H.Pagels and E. Tomboulis, Nucl. Phys. B143 (1978) 485; S.L.Adler, Phys. 
Rev. D33 (1981) 2905, D24 (1981) 1063(E), Phys. Lett. HOB (1982) 302, Nucl. 
Phys. B217 (1983) 381, S.L. Adler and T. Piran, Rev. Mod. Phys. 56 (1984) 1. 

[17] L.G. Yaffe, Rev. Mod. Phys. 541 (1982) 407, F.R. Brown and L.G. Yaffe, Nucl. 
Phys. B271 (1986) 267, T.A. Dickens et al, Nucl. Phys. B309 (1988) 1. 

[18] M. Baker, James S. BaU and F. Zachariasen, Phys. Rev. D37 (1988) 1036; D38 
(1988) 1926; D41 (1990) 2612. 

[19] C. Itzykson and J.B. Zuber, J. Math. Phys. 21 (1980) 411. 

[20] O. Haan, Phys. Lett. B106 (1981) 207; Yu. M. Makeenko, lectures in "Gauge 
Theories of the Eighties", Arctic School Proc, Acaslompolo, Finland, 1982, 
Eds. R. Ratio and J. Lindfors; P. Olesen, Nucl. Phys. B200 (1982) 381. 

[21] S. Levit, Strange Solutions of Mean Field Equations for Large N QCD in 1+1 Di- 
mensions, 28 International Winter Meeting of Nuclear Physics, Bormio (Italy), 
1990. 

[22] L.Salcedo, S.Levit and J.W.Negele, Nucl. Phys. B361 (1991) 585. 

[23] N.H. Crist and T.D. Lee, Phys. Rev. D22 (1980) 939; T.D. Lee, Particle Physics 
and Introduction to Field Theory, Harwood, New York, 1981, Ch.l8. 



28 



[24] R. Jackiw and J. Goldstone, Phys. Lett. 74B (1978) 81. 

[25] Yu.Simonov, Sov. J. Nucl. Phys. 41 (1985) 835,1014, Nucl. Phys. B307 (1988) 
512. 

[26] H.D. Dosch, Phys. Lett. 198B (1987) 177; H.D. Dosch and Yu. Simonov, Phys. 
Lett. 205B (1988) 339. 

[27] S. Mandelstam, Phys. Lett. 53B (1975) 476, Phys. Rep. 23C (1976) 245; t'Hooft, 
CERN report and papers quoted above, J. Kogut and L. Susskind, Phys. Rev. 
Dll (1975) 395. 

[28] LA. Batahn, S.G. Matinyan and G.K. Savvidi, Sov. J. Nucl. Phys. 26 (1977) 
214; Matinyan and G.K. Savvidi, B156 (1979) 1; G.K.Savvidy, Phys.Lett. 71B 
(1977) 133; 130B (1983) 303; Nucl. Phys. B246 (1984) 302; H. Pagels and E. 
Tomboulis, Nucl. Phys. B143 (1978) 485; B.Simon, Ann. Phys. 146, 209 (1983) 

[29] J.D.Bjorkcn, Elements of QCD; in: Lectures on Lepton-Nucleon Scattering and 
QCD; Prog. Nucl. Phys. Vol.4; eds.: W.B.Arwood, J.D.Bjorken, S.J.Brodsky 
and R.Straynowski. 

[30] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, Ch.3- 
6, McGrow-Hill,1965. 

[31] R. Bahan, Nuovo Cimento, B57 (1968) 183, and in Prom Microphysics To 
Macrophysics , vol. 1 , Springer - Verlag, Berlin, Heidelberg, 1991. 

[32] M.L.Mehta, Random Matrices, 2nd edition. Academic Press, 1991. 

[33] P. Ring and P. Schuck, The Nuclear Many-Body Problem, Springer- Verlag, 
New York, Heidelberg, Berlin, 1980. 



29 



